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I. INTRODUCTION 

In their seminal 1935 paper, Einstein, Podolsky, and 
Rosen (EPR) advocated that if "local realism" is taken 
for granted, then quantum theory is an incomplete de- 
scription of the physical world 1]. The EPR argument 
gained a renewed attention in 1964, when John Bell de- 
rived his famous inequalities, which must be satisfied 
within the framework of any local realistic theory 
The violation of Bell inequalities, predicted by quantum 
mechanics, has since then been observed in many exper- 
iments Illl[llEllllll[l3, thereby disproving the 
concept of local realism. So far, however, all these tests 
suffered from either a detector-efficiency loophole or a lo- 
cality loophole 111113, that is, the measured correlations 
may be explained in terms of local realistic theories ex- 
ploiting the low detector efficiency or the timelike interval 
between the two detection events 

A test of Bell inequality violation typically involves 
two distant parties, Alice and Bob, who simultaneously 
carry out measurements on parts of a shared quantum 
system that is prepared in an entangled state. Both par- 
ties randomly and independently decide between one of 
two possible quantum measurements oi, 02 and 61, 62. To 
avoid the locality loophole, the measurement events (in- 
cluding the choice of the measurement) at Alice's and 
Bob's sites must be spacelike separated. This suggests 
that optical systems are particularly suitable candidates 
for the test of Bell inequality violations. The technol- 
ogy of generation of entangled states of photons is very 
well mastered today and the prepared entangled states 
can be distributed over long distances via low-loss optical 
fibers (8|. However, the currently available single-photon 
detectors suffer from a too low efficiency 77, which opens 
the so-called detector-efficiency loophole. This loophole 
has been closed in a recent experiment with two trapped 
ions 0. However, the ions were held in a single trap, 
only several micrometers apart, so that the measurement 
events were not space-like separated. It was suggested 
that two distant trapped ions can be entangled via entan- 
glement swapping by first preparing an entangled state 



of an ion and a photon on each side and then projecting 
the two photons on a maximally entangled singlet state 
[TgHT^ lla . IT9I I . This technique could be used to close the 
locality loophole in the Bell test with trapped ions 1^. 
Very recently, the first step toward this goal, namely the 
entanglement between a trapped ion and a photon emit- 
ted by the ion, has been observed experimentally [20l |. 
However, the entanglement swapping would require in- 
terference of two photons emitted by two different ions, 
which is experimentally very challenging. 

An interesting alternative to the atom-based ap- 
proaches 2lL 22^ is represented by all-optical schemes 
involving balanced homodyne detection, which can ex- 
hibit very high detection efficiency |2j 01 • Unfortu- 
nately, the entangled two-mode squeezed state that can 
easily be generated experimentally |23| cannot 

be directly employed to test Bell inequalities with ho- 
modyning. This state is described by a positive defi- 
nite Gaussian Wigner function, which thus provides a 
local hidden variable model that can explain all corre- 
lations established via quadrature measurements carried 
by balanced homodyne detectors. Similarly as in the 
case of pur ification of continuous variable entanglement 
in ES in m HI, one has to go beyond the class of 
Gaussian states and Gaussian operations. For instance, 
it is possible to obtain a violation of Bell inequality with 
Gaussian two-mode squeezed vacuum state by perform- 
ing photon-counting measurements 'SS*] or the rather ab- 
stract measurements described in Refs. HE l36j . 
However, in contrast to balanced homodyning, these 
measurements are either experimentally infeasible or suf- 
fer from a very low detection efficiency. 

In order to close the detection loophole by using ho- 
modyne detectors, it is necessary to employ highly non- 
classical non-Gaussian entangled state whose Wigner 
function is not positive definite. Several recent theo- 
retical works indeed demonstrated that violation of Bell 
inequalities can be observed using balanced homodyn- 
ing m HE QEl' if specific entangled light states 
such as pair-coherent states, squeezed Schrodinger cat- 
like states, or specifically tailored finite superpositions of 
Fock states, are available. However, no feasible experi- 
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mental scheme is known that could generate the states 
required in Refs. jHIMIllliS- 

Recently, we have shown that a very simple non- 
Gaussian state obtained from two-mode squeezed vac- 
uumby subtracting a single photon fr'om each mode 
[43, can exhibit Bell violation with homodyning 
|4j|. An essential feature of our proposal is that the 
photon subtraction can be successfully performed with 
low-efficiency single-photon detectors, which renders the 
setup experimentally feasible. In fact, the basic building 
block of the scheme, namely the de-gaussification of a 
single-mode squeezed vacuum via single-photon subtrac- 
tion, has been recently successfully implemented experi- 
mentally jisj . 

In the present paper, we provide a thorough analysis of 
the scheme proposed in Ref. |3| • We present the details 
of the calculation of the Bell factor for a realistic setup 
that takes into account mixed input states, losses, added 
noise and imperfect detectors. Moreover, we shall also 
discuss several alternative schemes that involve the sub- 
traction of one, two, three, or four photons. The present 
paper is organized as follows. In Section II, we describe 
the proposed experimental setup and we introduce the 
Bell-CHSH inequalities. We then provide a simple pure- 
state analysis of the scheme assuming ideal detectors, 
which gives an upper bound on the achievable Bell viola- 
tion. In Section III, we present the mathematical descrip- 
tion of a realistic setup with imperfect detectors, losses 
and noise. Besides the scheme where a single photon is 
subtracted on each side, we will also analyze a scheme 
where two photons are subtracted on each side. This lat- 
ter scheme yields slightly higher Bell violation but only 
at the expense of a very low probability of state prepara- 
tion. Several other schemes composed of squeezed state 
sources, linear optics, and photon subtraction, are dis- 
cussed in Section IV. Finally, the conclusions are drawn 
in Section V. 



II. FEASIBLE BELL TEST WITH HOMODYNE 
DETECTION 

A. Proposed optical setup 

The conceptual scheme of the proposed experimental 
setup is depicted in Fig. 1. A source generates a two- 
mode squeezed vacuum state in modes A and B. This 
can be accomplished, e.g., by means of non-degenerate 
parametric amplification in a x^^^ nonlinear medium or 
by generating two single-mode squeezed vacuum states 
and combining them on a balanced beam splitter. Sub- 
sequently, the state is de-gaussified by conditionally sub- 
tracting a single photon from each beam. A tiny part of 
each beam is reflected from a beam splitter BSa (BS_b) 
with a high transmittance T. The reflected portions of 
the beams impinge on single-photon detectors such as 
avalanche photodiodes. A successful photon number sub- 
traction is heralded by a click of each photodetector PD^i 
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FIG. 1: Conceptual scheme of the proposed experimental 
setup for observing violation of Bell inequalities with balanced 
homodyning. The source emits two-mode squeezed vacuum 
mode in modes A and B. A small part of the beams is sub- 
tracted on two unbalanced beam splitters BS,4 and BSs and 
sent on single- photon detectors PDa and PDb- 



and PDb |43|. In practice, the photodetectors exhibit 
a single-photon sensitivity but not a single photon res- 
olution, that is, they can distinguish the absence and 
presence of photons but cannot measure the number of 
photons in the mode. Nevertheless, this is not a problem 
here because in the limit of high T, the most probable 
event leading to the click of a photodetector is precisely 
that a single photon has been reflected from the squeezed 
beam on the beam splitter. The probability of an event 
where two or more photons are subtracted from a single 
mode is smaller by a factor of « 1 — T and becomes to- 
tally negligible in the limit of T — » 1. Another important 
feature of the scheme is that the detector efficiency rj can 
be quite low because small r] only reduces the success rate 
of the conditional single-photon subtraction but it does 
not significantly decrease the fidelity of this operation. 
These issues will be discussed in detail in Section III. 

After generation of the non-Gaussian state, the two 
beams A and B together with the appropriate local os- 
cillators LOa and LOb are sent to Alice and Bob, who 
then randomly and independently measure one of two 
quadratures xf, , x^^ characterized by the relative phases 
61 , 02 and 01 , 02 between the measured beam and the 
corresponding local oscillator. The rotated quadratures 
= cos 9x^ + sin 9 and x^ = cos (j)x^ +sm(j) are 
defined in terms of the four quadrature components of 
modes A and B that satisfy the canonical commutation 
relations [x^ ,p^] = iSjk, j, k £ {A, B}. 

To avoid the locality loophole, the whole experiment 
has to be carried out in the pulsed regime and a proper 
timing is necessary. In particular, the measurement 
events on Alice's and Bob's sides (including the choice 
of phases) have to be space-like separated. A specific 
feature of the proposed setup is that the non-Gaussian 
entangled state needed in the Bell test is generated con- 
ditionally when both "event-ready" detectors 46] PD^i 
and PDb click. However, we would like to stress that 
this does not represent any loophole if proper timing is 
satisfied. Namely, in each experimental run, the detec- 
tion of the clicks (or no-clicks) of photodetectors PD^ 
and PD b at the source should be be space- like separated 
from Alice's and Bob's measurements. This guarantees 
that the choice of the measurement basis on Alice's and 
Bob's sides cannot in any way influence the conditioning 
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measurement or vice versa [1^ |4j, |46| ) . 

To demonstrate that the experimental data recorded 
by AUce and Bob are incompatible with the concept of 
local realism, we shall consider the Bell-CHSH inequal- 
ity originally devised for two-qubit system |47j . In this 
scenario, Alice (Bob) randomly and independently de- 
cides between one of two possible quantum measurements 
0-1, 0-2 (^1,^2) which should have only two possible out- 
comes +1 or —1. We define the Bell parameter 



S = (aibi) + (0162) + (0261) - (0262), 



(1) 



where (ajbk) denotes the average over the subset of ex- 
perimental data where Alice measured aj and, simultane- 
ously. Bob measured bk- If the observed correlations can 
be explained within the framework of the local-hidden 
variable theories, then S must satisfy the Bell-CHSH in- 
equality IS"! < 2. 

In the proposed experiment, Alice and Bob measure 
quadratures which have continuous spectrum. We dis- 
cretize the quadratures by postulating that the outcome 
is -1-1 when a; > and —1 otherwise. The two different 
measurements on each side correspond to the choices of 
two relative phases 0i, 6*2 and (j)i,4>2- Quantum mechani- 
cally, the correlation E{9j,(t)k) = {ajbk) can be expressed 
as 

/oc 
sign(a;^^ )P{xf^ , )dxj^ dx^ , (2) 

where P{x^^ ^^4,,) = i^e, > ^4,^ \Pc,ab\x^^ , x^, ) is the joint 
probability distribution of the two commuting quadra- 
tures xf, and x^^, and Pc^ab denotes the (normalized) 
conditionally generated non-Gaussian state of modes A 
and B. In practice, the correlations would be determined 
from the subset of the experimental data corresponding 
to the successful conditional de-gaussification, i.e., Alice 
and Bob would discard all results obtained in measure- 
ment runs where either PDa or PD^ did not click. We 
emphasize again that this does not open any loophole in 
the Bell test. 



B. Ideal photodetectors 

We shall first present a simplified description of the 
setup, assuming ideal photodetectors (??pd — 1) with 
single-photon resolution and conditionin g on detecting 
exactly a single photon at each detector |4ll |4^. This 
idealized treatment is valuable since it provides an upper 
bound on the practically achievable Bell factor S. More- 
over, as noted above, in the limit of high transmittance of 
BSa and BSb, T ^ 1, the reahstic (inefficient) detector 
with single-photon sensitivity is in our case practically 
equivalent to these idealized detectors. 

The two-mode squeezed vacuum state can be expressed 
in the Fock state basis as follows. 



|V'in(A))AS - Vl-A2^A"|n,n)AB, 



(3) 
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FIG. 2: (a) Bell factor 5* is plotted as a function of the effec- 
tive squeezing parameter TA for Oi = 0, O2 = 7r/2, 01 — — 7r/4 
and (j}2 = 7r/4. (b) Probability P of successful conditional 
generation of the state li/iout) as a function of the effective 
squeezing parameter TA, assuming T = 0.95. 



where A — tanh(s) and s is the squeezing constant. In the 
case of ideal photodetectors, the photon number subtrac- 
tion results in the state l-^out) aAflslV'inlrA)), where 
aA,B are annihilation operators and the parameter A is 
replaced by TA in order to take into account the trans- 
mittance of BSa and BSb- A detailed calculation yields 



out M-B 



(1-T2A^)3 
1 + T2A2 



^(n + l)(TA)"|n,n) AB, 



(4) 

and the probability of the conditional preparation of state 
Q can be expressed as 



1 + T^A" 
(1 -T2A2) 



3 • 



(5) 



For pure states exhibiting perfect photon-number cor- 
relations, the correlation coefficient depends only on 
the sum of the angles, E{9j,(j)k) = S{Oj -\- 4>k)- With 
the help of the general formula derived by Munro [s^ we 
obtain for the state 



f(^) 



(1 -T2A2)3 



E 



87r(2TA) 



1 + T^A^ n\m\(n — m)'^ 

n>rn ^ ' 



(n -I- l)(m -I- 1) 



X \T(n^ m) — !F{m, n)]^ cos[(n — m)(y9]. 



(6) 



n=0 



where T{n,m) = T-^{{1 - n)/2)T-^{-m/2) and T{x) 
stands for the Euler gamma function. 

We have numerically optimized the angles 9i^2 and 0i^2 
to maximize the Bell factor S. It turns out that for any 
A, it is optimal to choose 61 = 0, 82 = 7r/2, = —tt/A 
and (j)2 — tt/4:. The Bell factor S for this optimal choice 
of angles is plotted as a function of the effective parame- 
ter TA in Fig. 2(a), and the corresponding probability of 
success of the conditional preparation of the state IV'out) 
is plotted in Fig. 2(b). We can see that S is higher 
than 2 so the Bell inequality is violated when TA > 0.45. 
The maximal violation is achieved for TA ~ 0.57, giving 
S ~ 2.048. This figure is quite close to the maximum Bell 
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factor S = 2.076 that could be reached with homodyne 
detection, sign binning, and arbitrary states exhibiting 
perfect photon-number correlations \ip) = J2n'^n\n,n) 



III. REALISTIC MODEL 

In this section we will consider a realistic scheme 
with inefficient (rypD < 1) photodetectors exhibiting sin- 
gle photon sensitivity but no single-photon resolution, 
and realistic homodyning with efficiency ryBHD < 1- 
The mathematical description of this realistic model of 
the proposed experiment becomes strikingly simple if 
we work in the phase-space representation and use the 
Wigner function formalism. Even though the state used 
to test Bell inequalities is non-Gaussian, it can be ex- 
pressed as a linear combination of four Gaussian states, 
so all the powerful Gaussian tools can still be used. 

This section is further divided into three sub-sections. 
The first one gives a brief overview of the Gaussian states, 
linear canonical transformations of quadrature operators, 
and Gaussian completely positive maps. In the second 
subsection, an analytical formula for the Bell factor S 
is derived, and the influence of detector inefhciencies, 
losses, and noise on the proposed Bell experiment is in- 
vestigated in detail. Finally, an extended setup involving 
two-photon subtraction from each mode is studied in the 
third subsection. 



A. Gaussian states and Gaussian operations 

In quantum optics Gaussian states are often encoun- 
tered as states of n modes of light. These states are com- 
pletely specified by the first and second moments of the 

quadrature operators with r = {xi,pi, ,Xn,Pn)^- 

Here satisfy the canonical commutation relations 
(CCR) [xj,pk] = iSj^k- Instead of referring to the density 
matrix one may refer to the Wigner function defined on 
phase space 



W 



1 



7r"-\/det 7 



exp[ 



^{r-dfj-\r-d)], (7) 



where d is the vector of first moments, dj — (rj), and 7 
is the covariance matrix 



= {nrj + rjn) - 2didj 



(8) 



In this paper we shall deal only with states with zero dis- 
placement dj = 0. Some relevant examples of Gaussian 
states that we shall need in what follows include: (i) The 
n-mode vacuum state with dj = and covariance matrix 
equal to the identity matrix, 7vac — hn ■ (ii) Single- mode 
squeezed vacuum state with dj — and covariance ma- 
trix 



7SMS 



=2s 





-2s 



(9) 



where s is the squeezing parameter, (iii) Two-mode 
squeezed vacuum state with dj — and covariance ma- 
trix 



7TMS 



cosh(2s) 


sinh(2s) 






cosh(2s) 


— sinh(2s) 



sinh(2s) 


cosh(2s) 






— sinh(2s) 


cosh(2s) 



(10) 

Optical operations that can be implemented with beam 
splitters, phase shifters, squeezers and homodyne detec- 
tion correspond to Gaussian operations. Their important 
property is that they map a Gaussian input state onto 
a Gaussian output state. Gaussian unitary transforma- 
tions realize the mapping r r' = Sr which preserves 
the CCR. This is the case if S' € Sp{2n, R) the so-called 
real symplectic group. On the covariance matrix level 
the transformation reads 



7 ^ SjS^ 



(11) 



A particular subset of symplectic transformations is 
formed by the symplectic matrices S that are also or- 
thogonal S G Sp{2n, E.) n 0(2n). Those transformations 
are called passive because they do not change the total 
number of photons. The most common passive trans- 
formations include (i) mixing two modes of light with 
a beam splitter of (intensity) transmittance T and re- 
flectance 1 — T 



'BS 



Vt 








T 















-x/i~T 



T 
Vt 



and a phase shift of a single mode 



SpsiO) = 



cos 6 sin t 
— sin 9 cos ( 



(12) 



(13) 



All passive linear canonical transformations of n modes 
can be implemented by optical interferometers consisting 
of beam splitters and phase shifters. 

The second group of linear canonical transformations 
are the active transformations that describe phase sensi- 
tive amplification of light. The archetypal examples are 
the single-mode squeezer 



SsMS = 



and the two-mode squeezer 





e-' 



(14) 



'S'tms — 
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sinh(s) 






cosh(s) 


— sinh(s) 



sinh(s) 


cosh(s) 






— sinh(s) 


cosh(s) 



(15) 
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These matrices describe the operation of ideal degenerate 
(5sMs) or nondegenerate (^tms) optical parametric am- 
plifier (OPA). In particular, a nondegenerate OPA pro- 
vides a source of entanglement since it transforms the 
input vacuum into a two-mode squeezed vacuum state. 

Noisy channels and phase-insensitive amplifiers are ir- 
reversible quantum operations which cannot be described 
by Gaussian unitary transformations. Instead, they 
can be modeled within the more general framework of 
trace-preserving Gaussian completely positive (CP) maps 
1281 129II . The covariance matrix transformation reads 



(ideal) detectors is related to via a Gaussian CP map. 



7 ^ A-fA^ + G. 



(16) 



Of particular importance is the propagation through a 
lossy quantum channel with transmittance 77, which is 
characterized hy A = ^ I and G = (1 — 77)/. In what 
follows, we shall use lossy channels followed by perfect 
detectors to model inefficient detectors. 



B. Two photon subtractions 

We shall now present a detailed calculation of the Bell 
factor for the setup depicted in Fig. 1. Our model takes 
into account realistic photodetectors (yypD < 1) with 
single-photon sensitivity, imperfect homodyning and the 
added electronics noise. 



1. Preparation of a non- Gaussian state 

As shown in Fig. 1, the modes A and B are ini- 
tially prepared in a two-mode squeezed vacuum state, 
and the auxiliary modes C and D are in vacuum state. 
The Wigner function of the four-mode state ABCD is a 
Gaussian centered at the origin. 



1 



in, ABCD 



7r^-\/det7i: 



: exp [-r^7i„V] 



(17) 



where r — [x^jp^^, . . . , ,p^]. The initial state is fully 
characterized by the covariance matrix 



where 



S,^ — y/7]BiiD I AB ' 



G = (1 - r]BBB)lAB © (1 - ■'Jpb)Icd, 
and the symplectic matrix 



BS,AC ' 



BSMD 



(19) 



(20) 
(21) 



(22) 



describes the mixing of modes A with C and B with D 
on the unbalanced beam splitters BSa and BSb, respec- 
tively. 

The state Pc,AB is prepared by conditioning on observ- 
ing clicks at both photodetectors PD^i and PD^. These 
detectors respond with tho different outcomes, either a 
click, or no click. Mathematically, an ideal detector 
with a single photon sensitivity is described by a two- 
component positive operator valued measure (POVM) 
consisting of the projectors onto the vacuum state and on 
the rest of the Hilbert space, no = |0)(0|,ni=/-|0)(0|. 
The resulting conditionally prepared state PcAB can be 
calculated from the density matrix Po\it,ABCD as follows, 

Pc,AB = ^T:CD[Pont,ABCD{lAB ® IIi^c '^1,d)]- (23) 

It is instructive to rewrite the partial trace in Eq. H23() 
in terms of Wigner functions, taking into account that 



/oo 
Wx{r)WY{r)S^r, 
-OO 



(24) 



7in = 7TMS,AB 



CD, 



(18) 



where Wxir) and IlV(r) denote the Wigner representa- 
tions of the operators X and Y , respectively, and N is 
the number of modes we trace over. The POVM element 
Hi is a difference of two operators whose Wigner rep- 
resentations are both Gaussian functions, Wi = l/(27r). 
Wo = 7r~^e~^ ~P . After a bit lengthy but otherwise 
straightforward calculations we find that the Wigner 
function Wc,ab of (normalized) conditionally prepared 
state (|23|) can be expressed as a linear combination of 
four Gaussian functions, 



where 7tms is the covariance matrix of a two-mode 
squeezed vacuum (fTT)| and denotes the direct sum of 
matrices. 

The imperfect single-photon detectors (balanced ho- 
modyne detectors) with detector efficiency ?7pd ('7bhd) 
are modeled as a sequence of a lossy channel with trans- 
mittance ?7pd (?7bhd) followed by an ideal photodetector 
(homodyne detector). In our setup, the modes AC (BD) 
interfere on the unbalanced beam splitters BSa (BSb) 
and pass through the four "virtual" lossy channels before 
impinging on ideal detectors. The covariance matrix of 
the mixed Gaussian state Pout.ABCD just in front of the 



Wc,AB{r) = 



TT^ Pg VdetTout ^ ^det Tj^cd 



^-r'^T^_ABr 



(25) 

where qi = 1, (72 = 93 = —2 and 54 = 4. The correspond- 
ing probability of success is given by 



Pa 



1 



Vdet 7out ^ ^/A.ei{Tj^ABTj,cD) 



(26) 



To define the various matrices appearing in Eqs. H25|l 
and (|26|) . we first introduce a matrix F = 7^^ and we 
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divide T into four smaller submatrices with respect to 
the bipartite AB vs CD splitting, 



r = 



Tab cr 

r 



CD 



It holds that 

^j,AB — Tab 
and the four matrices Tj cd read 



^i.CD = rcn, 

^2,CD — ^CD + Ic ® Od, 

^s.CD — TcD + 0c ® Id, 
^4,CD = Tcu + /cr>- 



(27) 



(28) 



(29) 



2. 



Correlation coefficient E{9j,(j)k) 

The joint probability distribution P{x-^,,x^^) of the 
quadratures xf. and x^^ appearing in the formula (|2J) for 
the correlation coefficient E{dj, (/)fc) can be obtained from 
the Wigner function as a marginal distribution. We 
have 



W, 



c.AB 



{Slre,,4>.)dpldpl, (30) 



— oo — oo 



] and the symplectic ma- 
fc) describes local phase 



where rg^ = [x^^ , pf^ , x^^ , p; 

trix Ssh = S'ps,A(6'i) © Sps^B 
shifts applied to modes A and B that map the measured 
quadratures xf, and x^^ onto the quadratures x^ and 

x^ , respectively. 

In order to express the result of the integration in Eq. 
()30|l in a compact matrix notation, we re-order the ele- 
ments of the vector rg.^^^, as follows, 



(31) 



which defines a matrix 5'hom- After these algebraic ma- 
nipulations, the four matrices ^j,AB appearing in the ex- 
ponents in Eq. I|25|l transform to 
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^j,AB — ShomSsh^ j,ABS'[i^S^^^^-^ — 



^0 Q 



(32) 



where we have divided the matrix into four sub- 

matrices with respect to the a; vs p splitting. A straight- 
forward integration over pf, and p^^ in Eq. 130(1 then 
yields the joint probability distribution. 



1 



7rPG\/det7out ^ detV j d \/ Act B j ' 

(33) 



where y = {x'^.^x^^Y' ^-i^d 

T,=A,-C,B-'CJ. 



(34) 



Taking into account the choice of binning, the normal- 
ization of the joint probability distribution, and its sym- 
metry, P{xg^_,x? ) = P{—xf,, —xf ), we can express the 
correlation coefficient as follows, 

E{e,,^u) = 4 / / P{x^^,xl)dx^^dxl - 1. (35) 
Jo Jo 

This last integral can be easily evaluated analytically. 
For a given V j matrix 



(36) 



the integral of the exponential term 



Gj = / / e~''^^i~^^^2-2c,«iy2^y^rfy2 (37) 
Jo Jo 

can be calculated by transforming to polar coordinates 
and integrating first over the radial coordinate and then 
over the angle. After some algebra, we finally arrive at 



1 



2i / a,-5, — c 



arctan 

2 



(38) 



The final fully analytical formula for the correlation co- 
efficient reads 



7rPG\/det7out 



E 



\ ^/det Vj.cD ^/detSj 



-1 



(39) 



and the Bell factor can be expressed as 
S = S(0i, (/-i) + E(d^,<^2) + E(B^, ^i) - E(Q^, 4>2). (40) 

3. Violation of Bell-CHSH Inequalities 

A necessary condition for the observation of a viola- 
tion of Bell inequalities with homodyne detectors is that 
the Wigner function of the two-mode state used in the 
Bell test is not positive definite. Figure illustrates that 
the Wigner function ((25(1 of the conditionally generated 
state Pc,AB is indeed negative in some regions of the phase 
space. The area of negativity, as well as the attained neg- 
ative values of W, are rather small, which indicates that 
we should not expect a high Bell violation with homo- 
dyning. 

As we have shown in Section II, the maximum Bell fac- 
tor S achievable with our setup and sign binning is about 
5 = 2.048. We conjecture that this binning is optimal or 
close to optimal. This is supported by the simple struc- 
ture of the joint probability distribution 133() . As can be 
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seen in Fig.2fa,b), P exhibits two peaks, both located in 
the quadrants where AHce's and Bob's measured quadra- 
tures have the same sign. Note also that the two-peak 
structure is a clear signature of the non-Gaussian char- 
acter of the state (c.f. Fig. 2fc,d)). We have carried out 
numerical calculations of S for several other possible bin- 
nings which divide the quadrature axis into three or four 
intervals, and have not found any binning which would 
provide higher S than the sign binning. We have also 
performed optimization over the angles 9j and ipk and all 
the results and figures presented in this Section were ob- 
tained for the optimal choice of angles 6i = 0, 6*2 = 7r/2, 
01 = -7r/4, (j)2 = 7r/4. 

Figure Ufa) illustrates that the Bell-CHSH inequality 
IS*! < 2 can be violated with the proposed set-up, and 
shows that there is an optimal squeezing Aopt which max- 
imizes S. This optimal squeezing is well predicted by 
the simple model assuming perfect detectors with single- 
photon resolution (section II B), Aopt?" ~ 0.57. The curve 
plotted for T — 0.99 practically coincides with the results 
obtained from the simple model presented in Sec. II B, 
c.f. Fig. 2(a). This confirms that in the limit T ^ 1 
the detectors with single-photon sensitivity become for 
our purposes equivalent to photodetectors with single- 
photon resolution. The maximum Bell factor achievable 
with our scheme is about S'max ~ 2.045 which represents 
a violation of the Bell inequality by 2.2%. To get close to 
the S'max one needs sufficiently high (but not too strong) 
squeezing. In particular, the value A ~ 0.57 corresponds 
to approximately 5.6 dB of squeezing. Figure |Sfb) illus- 
trates that there is a clear trade-off between S and the 
probability of success Pq ■ To maximize 5* one should use 
highly transmitting beam splitters but this would reduce 
Pg- The optimal T that should be chosen would clearly 
depend on the details of the experimental implementa- 
tion. 

It follows from Fig. Ela) that the Bell factor S de- 
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FIG. 3: A one-dimensional cut of the Wigner function of the 
two- mode state pc,AB along the line = x^, — = 
for A = 0.6 and beam splitters BSa and BSb transmittances 
T = 0.95. Notice the regions where W is negative. 




FIG. 4: Joint probability distribution P{x^. , x^^). Panels (a) 
and (b) show the distribution for the conditionally-prepared 
non-Gaussian state with T = 0.99. Panels (c) and (d) display 
the distribution for the initial Gaussian two-mode squeezed 
vacuum state. The curves are plotted for perfect detectors 
ripo — rjBHD = 100%, squeezing A = 0.6 and ^auco = and 
0Bob = 7r/4. 
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FIG. 5: Violation of Bell-CHSH inequality with the 
conditionally-prepared non-Gaussian state, (a) Bell factor 5* 
as a function of the squeezing, (b) Probability of success 
of the generation of the non-Gaussian state as a function of 
the squeezing. The curves are plotted for perfect detectors 
(??PD = riBHD = 100%) with T = 0.9 (solid line), T = 0.95 
(dashed line), and T — 0.99 (dot-dashed line). 

pends only very weakly on the efficiency rypD of the single- 
photon detectors, so the Bell inequality can be violated 
even if yypD ~ 1%. This is very important from the ex- 
perimental point of view because, although the quantum 
detection efficiencies of the avalanche photodiodes may 
be of the order of 50%, the necessary spectral and spatial 
filtering which selects the mode that is detected by the 
photodetector may reduce the overall detection efficiency 
to a few percent. Low detection efficiency only decreases 
the probability of conditional generation Pq of the non- 
Gaussian state, see Fig.EJb). The dependence of Pq on 
77PD and T can be very well approximated by a quadratic 
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FIG. 6: Effect of the inefficiency of tfie pfiotodetectors PDa 
and PT>B ■ (a) Bell parameter S as a function of the efficiency 
rjpu of the photodetectors. (b) Probability of success as a 
function of the efficiency r;pD- The curves are plotted for 
T\ = 0.57, r/BHD ~ 100% and the same transmittances as in 
Fig. El 
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FIG. 7: Effect of inefficient homodyning. (a) Bell parameter 
5 as a function of the efficiency 77BHD of the homodyning. 
The curve is plotted for TA = 0.57, r]po = 30% and the same 
transmittances as in Fig. |3 (b) Bell parameter achieved for 
the optimal squeezing Aopt is plotted as a function of r;BHD. 
(c) Optimal squeezing Aopt is plotted as a function of r?BHD. 
The curve is plotted for rjpo = 30% and the same transmit- 
tances as in Fig. |3 

function, « Vpui^ — T)^ which quickly drops when 
?7pd decreases. In practice, the minimum necessary rjpD 
will be determined mainly by the constraints on the to- 
tal time of the experiment and by the dark counts of the 
detectors. 

In contrast, the Bell factor S strongly depends on the 
efficiency of the homodyne detectors, and ?7bhd must be 
above ~ 90% in order to observe Bell violation, see Fig. 
[7| However, this is not an obstacle because such (and 
even higher) efficiency has been already achieved exper- 
imentally (see e.g. |43|). Interestingly, we have found 
that it is possible to partially compensate for imperfect 
homodyning with efficiency ?7bhd < 1 by increasing the 
squeezing of the initial state. This effect is illustrated in 
Fig. [7Jb) which shows the dependence of the Bell factor 
S on ?7bhd for optimal squeezing Aopt. Figure EJc) then 
shows how the optimal squeezing increases with decreas- 
ing ?7bhd- 

In addition to imperfect detection efficency 77BHD, the 
electronic noise of the homodyne detector is another fac- 
tor that may reduce the observed Bell violation. We 
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FIG. 8: Effect of the electronic noise and thermal input states, 
(a) Maximum achievable Bell parameter S with the optimal 
squeezing Aopt as a function of the electronic noise A'oi. (b) 
Optimal squeezing Aopt giving the highest Bell parameter S 
for a given electronic noise, (c) Maximum Bell parameter S 
as a function of the thermal noise of the input state Vnoisc. 
(d) Optimal squeezing Aopt giving the highest Bell parameter 
S for a given thermal noise at the input. The curves are 
plotted for 77PD = 30%, 77BHD = 95%, and T = 0.9 (solid 
line), T = 0.95 (dashed line), and T = 0.99 (dot-dashed line). 

model the added electronic noise by assuming that the 
effective quadrature that is detected Xdct is related to the 
signal quadrature xs by a formula, 

^noisc , 

where x^ac and Xnoiso are two independent Gaussian dis- 
tributed quadratures with zero mean and variance 1/2, 
and TVoi is the electronic noise variance expressed in shot 
noise units. On the level of covariance matrices, A'oi can 
be included by modifying the formula for the noise matrix 
G, 

G = (1 - r?BHD + N^i)Iab © (1 - Vpb)Icd- (41) 

The homodyne detector with electronic noise is actually 
equivalent to a detector without noise but with a lower 
homodyne detector efficiency ti'^ylu ~ '?bhd/(1 + -^ci)- 
This can be shown by noting that the re-normalized 
quadrature Xdot / V 1 + -^oi is exactly a quadrature that 
would be detected by a balanced homodyne detector with 
A'ci = and efficiency ?7bhd ■ calculations reveal that 
the electronic noise should be 15 — 20 dB below shot noise 
(see Fig. |HIa) and (b)), which is currently attainable with 
low-noise charge amplifiers. Again, higher squeezing can 
partially compensate for the increasing noise. 

So far we have assumed that the source in Fig. 1 emits 
pure two-mode squeezed vacuum state. However, experi- 
mentally, it is very difficult to generate pure squeezed vac- 
uum saturating the Heisenberg inequality. It is more real- 
istic to consider a mixed Gaussian state such as squeezed 
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thermal state which can be equivalently represented by 
adding quadrature independent Gaussian noise with vari- 
ance Koiso to each mode of the two-mode squeezed vac- 
uum. The effect of the added noise stemming from input 
mixed Gaussian state is quite similar to the influence of 
the electronic noise of the homodyne detector, see Fig. |H| 
(c) and (d). We find again that the added noise in the 
initial Gaussian state should be 15 — 20 dB below the 
shot noise. 



C. Four photon subtractions 

Until now we have focused on a single-photon subtrac- 
tion on each side (one photon removed from mode A and 
one from mode B). If we now consider a scheme where 
two photons are subtracted from each mode, the de- 
gaussification of the state will be stronger and we may ex- 
pect a higher Bell violation than before. To subtract two 
photons from each mode, we only need to add one more 
unbalanced beam splitter and photodetector on each side 
in Fig. 1. A successful state generation would be indi- 
cated by simultaneous clicks of all four detectors. Assum- 
ing perfect photon-number resolving detectors, the state 
generated from two-mode squeezed vacuum (O by sub- 
tracting two photons from each mode can be expressed 
as 

\i>out) AB OC a^6^|l/'in(T^A))AS 

(X J2{n + 2){n+l){T^Xr\n,n)AB,m 

n 

and the probability of success reads 



P4ph = 2T2(l-r)4A'*(l-A2 



1 + lOT^X^ + T^X^ 
' (1-T4A2)5 ■ 



(43) 



Since the state (I42f) exhibits perfect photon number cor- 
relations, the Munro's formula for the Bell factor can 
again be directly applied 39]. Numerical calculations 
show that the maximum Bell violation with the state 
(|42|l and sign binning of quadratures is achieved for 
T^A = 0.40 which yields S'max,4ph = 2.064, which is 
indeed higher than the maximum achievable with two- 
photon subtraction, S'max,2ph = 2.048, and very close to 
the maximum value S — 2.076 |39j . 

A more realistic description of the four-photon sub- 
traction scheme that takes into account realistic imper- 
fect detectors can be developed using the approach de- 
scribed in detail in Sec. IIIB. We find that the Wigner 
function of the conditionally generated state is a linear 
combination of sixteen Gaussians. The results of numeri- 
cal calculations are shown in Figs. El^a) and (b), which il- 
lustrate that the two-photon subtraction from each mode 
yields higher violation of Bell-CHSH inequality than one- 
photon subtraction only for very high transmittances 
T > 0.95. For lower transmittances, the fact that the 
photodetectors do not distinguish the number of pho- 
tons reduces the Bell factor. Moreover, adding a second 
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FIG. 9: Violation of Bell-CHSH inequality with four pho- 
ton subtractions, (a) Bell parameter 5 as a function of the 
squeezing A for perfect detectors rypD ~ t/bhd ~ 100%. (b) 
Bell parameter S' as a function of the efficiency i]bho of the ho- 
modyning. The curve is plotted for T^A = 0.40, r/pn = 100% 
and the same transmittances as in Fig. |^ 



a) 



c) 



b) 



d) 



FIG. 10: Symbol convention, (a) Single-mode squeezer along 
the X quadrature, (b) Two-mode squeezer, (c) Beam splitter, 
(d) Conditional subtraction of a photon as described in the 
preceding section. 



stage of photon subtractions dramatically decreases the 
probability of generating the non-Gaussian state. The 
probability can be estimated as Pq ~ Vpui^ ^ T)"^, so 
for T > 0.95 and 77 = 50% we get Pq « 10"^ and the 
duration of data acquisition would make the experiment 
infeasible. We conclude that from the practical point of 
view there seems to be no advantage in using the scheme 
with four photon subtractions instead of the much sim- 
pler scheme with two photon subtractions. 



IV. ALTERNATIVE SCHEMES 

In this section we will study the violation of Bell-CHSH 
inequalities for a large group of alternative schemes, 
which involve from one to four photon subtractions. The 
main objective of this section is to compare the maximum 
Bell-CHSH factor S obtained for the different proposed 
setups. As the main purpose of this section is the com- 
parison of the different schemes, we will consider only 
idealized schemes with almost perfect single-photon sub- 
traction on the beam splitters (T = 0.99), and perfect 
photodetectors and homodyning (jypD = '7bhd — 100%). 
The maximum achievable Bell factor for each scheme pre- 
sented below was determined by optimizing over the an- 
gles 9i^2, 4'i,2 as well as over the squeezing A of the ini- 
tial Gaussian states. The sign binning of the measured 
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Schemes: one subtraction 
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FIG. 11: Schemes with only one photon subtraction. The first 
column labels the different setups proposed, the second shows 
the scheme and finally the last column gives the maximal Bell 
factor S. (a) Photon subtraction after the creation of the 
two-mode squeezed vacuum, (b) Photon subtraction before 
mixing two single-mode squeezed states on a beam splitter. 



quadratures has been used in all cases. All the schemes 
presented in this section use the symbol convention de- 
picted in Fig. ^1 

In the preceding section, we have seen that the prob- 
ability of successful generation of a non-Gaussian state 
decreases significantly with the number of photon sub- 
tractions. At the same time the complexity of the imple- 
mentation of the experimental setup increases with the 
number of photon subtractions. It is then obvious that 
the most interesting schemes for a Bell-CHSH violation 
are those involving only one photon subtraction. Unfor- 
tunately, for the schemes that we have considered (see 
Fig. Ill|) . no violation is observed. 

After one photon subtraction, the simplest schemes are 
those with two photon subtractions. In the preceding 
sections it was shown that it is possible to violate Bell- 
CHSH inequality with two photon subtractions (scheme 
Fig. I12r a')'). It follows from Fig. E|that several other 
schemes (scheme Fig. I12f d') and (e)) also violate Bell- 
CHSH inequality but the maximal achievable Bell factor 
S is much smaller in comparison to the scheme shown in 
Fig. da). 

By adding one more photon subtraction to the schemes 
shown in Fig. 1121 we can construct an ensemble of 
schemes with three photon subtractions. After numerical 
optimization we have found that none of these schemes 
succeeds to violate Bell-CHSH inequality. This striking 
result together with the the fact that we have not found 
any violation for schemes based on a single subtraction 
suggests that it may be necessary to have a scheme with 
an even number of photon subtractions in order to ob- 
serve S > 2. 

In the preceding section, we have also proposed one 
scheme with four photon subtractions that violates Bell- 
CHSH inequality. Many other possible schemes exist 
where four photons are subtracted. Figure [TSl illustrates 
some particular examples, which are based on the prepa- 
ration of two-mode squeezed vacuum via mixing of two 
single-mode squeezed states on an balanced beam split- 
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Schemes: two subtractions 
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FIG. 12: Schemes with two photon subtractions. The right 
column gives the maximal value of the Bell factor S for the 
proposed setups. 
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Schemes: four subtractions 
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FIG. 13: Schemes with four photon subtractions. Last col- 
umn gives the maximal value of the Bell factor 5 for the 
proposed setups. 



ter. The photon subtractions are symmetrically placed 
to both modes. Strikingly, if all four photons are sub- 
tracted either before or after mixing on a beam splitter, 
then we get S > 2. However, if a single photon is sub- 
tracted from each mode both before and after combining 
the modes on a beam splitter, then we do not obtain any 
Bell violation. 

Finally we have also studied an alternative group of 
schemes where instead of subtracting photons separately 
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from modes A and B, we mix the auxiliary modes C 
and D on a balanced beam splitter before the detection 
on the photodetectors. Consider the scheme depicted 
in Fig. I14f a) where only a single photon is subtracted. 
The mixing of modes C and D on a beam splitter erases 
the information about the origin of the detected photon 
which implies that the conditionally prepared state is a 
coherent superposition of states where a single photon 
has been removed either from mode A or from mode B. 
However, even this modification does not lead to Bell 
violation with just a single subtraction. 

We can extend the scheme by placing a photodetector 
at both output ports of the beam splitter, cf. Fig. IHT b). 
In the limit of a high transmittance T — > 1, the condi- 
tioning on the click of each detector selects the events 
where there were altogether two photons at the beam- 
splitter inputs. The bosonic properties of the photons 
imply that a simultaneous click of both photodetectors 
occurs only if the two subtracted photons are coming 
from the same mode (A or B) ji^, but again we do not 
know from which mode the two photons are subtracted. 
This scheme is thus equivalent to the superposition of 
two schemes of the type shown in Fig. Ifif c') . Unlike the 
scheme in Fig. I12r c'). the scheme in Fig. I14r b') is sym- 
metric with respect to the modes A and B. However, 
no violation can be observed. On the other hand, the 
scheme in Fig. ICT c) leads to S' > 2 by realizing a super- 
position of states where two photons are subtracted from 
a single-mode squeezed vacuum state and this state is 
then mixed with another single-mode squeezed vacuum 
on a balanced beam splitter, see ll2f d). In comparison to 
the scheme in Fig. I12r d). we obtain much higher viola- 
tion S = 2.046. 



V. CONCLUSIONS 

We have proposed an experimentally feasible setup al- 
lowing for a loophole-free Bell test with efficient homo- 
dyne detection using a non-gaussian entangled state gen- 
erated from a two-mode squeezed vacuum state by sub- 
tracting a single photon from each mode. We have pre- 
sented a full analytical description of a realistic setup 
with imperfect detectors, noise and mixed input states. 
We have studied in detail the influence of the detec- 
tor inefficiencies, the electronic noise of homodyne de- 
tector, and the input mixed states, on the achievable 
Bell violation. The main feature of the present scheme is 
that it is largely insensitive to the detection efficiency of 
the avalanche photodiodes that are used for conditional 
preparation of the non-gaussian state, so that detector 
efficiencies of the order of a few per cent are sufficient. 
On the other hand, the detection efficiency of the bal- 
anced homodyne detector should be of the order of 90% 
and the electronic noise of the homodyne detector should 
be at least 15 dB below the shot noise level. The optimal 
squeezing that yields maximum Bell violation depends on 
the experimental circumstances but is, generally speak- 



Schemes: Superposition S 




FIG. 14: Schemes consisting of superpositions of other 
schemes proposed above, (a) Superposition of one photon 
subtraction on modes A ot B. (b), (c) Superposition of two 
photon subtractions on modes A or B. 

ing, within the range of experimentally attainable values. 
As a rule, the optimal squeezing increases with decreas- 
ing ryBHD and increasing noise. 

We have also discussed several alternative schemes that 
involve the subtraction of one, two, three or four photons. 
The experimentally simplest and most appealing schemes 
are those where only a single photon is subtracted be- 
cause photon subtraction is a delicate operation and also 
each subtraction in the scheme drastically reduces the 
probability of successful state generation. Unfortunately, 
we have not been able to find a scheme with only a sin- 
gle subtraction which would exhibit violation of Bell in- 
equalities. However, the class of schemes that we have 
studied is still somewhat restricted. One can thus hope 
that such a scheme may be designed by considering more 
complicated setups involving unbalanced beam splitters 
and possibly a different binning procedure |53- This is- 
sue certainly deserves further investigation. 

Among all the schemes where two photons are sub- 
tracted, the maximum violation S = 2.046 is achieved by 
the scheme discussed in Sections H and HI. Taking into 
account that we have not found any scheme with three 
photon subtractions which would violate BcU-CHSH in- 
equality, the only way of exceeding the 2.046 violation 
appears to be by subtracting four photons. This scheme 
has been analyzed in some detail in Sec. IIIC where it 
was shown that this allows us to reach the Bell factor 
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S = 2.06. Unfortunately, the price to pay for this shght 
increase of S is that the probabihty of successful condi- 
tional generation is so low that it makes the experiment 
infeasible. 

The results presented in this paper provide a clear 
example of the utility of conditional photon subtrac- 
tion which can be considered as an important novel tool 
in quantum optics and quantum information processing 
with continuous variables. Besides violation of Bell in- 
equalities, this method can be used to generate highly 
non-classical states of light to improve the fide lity 
of teleportation of continuous variable states |4lL 142. |43 | 
and it forms a key ingredient of recently proposed entan- 
glement purification protocols for continuous variables 
|3lL Is^ . The very recent experimental demonstration of 
a single photon subtraction from a single-mode squeezed 
vacuum state provides a strong incentive for further theo- 
retical and experimental developments along these lines, 
and we can thus expect that some of the schemes dis- 
cussed in the present paper will be experimentally im- 



plemented in a not too distant future. 

Note added: After this work was completed we have 
learned that a scheme for observing a violation of Bell 
inequalities similar to the scheme discussed in Sec. II of 
the present paper has been independently proposed by 
H. Nha and H.J. Carmichael HJ. 
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